Anyons in a weakly interacting system 
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In quantum theory, indistinguishable particles in three-dimensional space behave in only two 
distinct ways. Upon interchange, their wavefunction maps either to itself if they are bosons, or 
to minus itself if they are fermions. In two dimensions a more exotic possibility arises: upon 
exchange of two particles called "anyons" the wave function acquires phase e'° 7^ ±1. Such fractional 
exchange statistics are normally regarded as the hallmark of strong correlations. Here we describe 
a theoretical proposal for a system whose excitations are anyons with the exchange phase 9 — 7r/4 
and charge — e/2, but, remarkably, can be built by filling a set of single-particle states of essentially 
noninteracting electrons. The system consists of an artificially structured type-II superconducting 
film adjacent to a 2D electron gas in the integer quantum Hall regime with unit filling fraction. 
The proposal rests on the observation that a vacancy in an otherwise periodic vortex lattice in 
the superconductor creates a bound state in the 2DEG with total charge —e/2. A composite of 
this fractionally charged hole and the missing flux due to the vacancy behaves as an anyon. The 
proposed setup allows for manipulation of these anyons and could prove useful in various schemes 
for fault-tolerant topological quantum computation. 



I. INTRODUCTION 

Anyons and fractional charges appear in a variety of 
theoretical models involving electron or spin degrees of 
freedom in 2D [J H i, 0, H @. In all known cases, in- 
teractions between the fundamental degrees of freedom 
are of paramount importance for the emergence of the 
exotic excitations. For example, it is understood that 
the Coulomb interaction between electrons is crucial to 
stabilize the incompressible ground state in the frac- 
tional quantum Hall (FQH) liquids that support frac- 
tionalized excitations. One can say that such systems 
are strongly correlated in the sense that their many-body 
wavcfunctions cannot be constructed from the single- 
particle states of the constituent fundamental degrees of 
freedom. 

The question thus arises whether anyons and frac- 
tionalization are inextricably connected with the phe- 
nomenon of strong correlations as defined above. There 
are compelling reasons to believe that this is in fact 
not so. In one dimension, Su, Schrieffer and Heeger [7( 
showed that a half of an electron can be bound to a do- 
main wall in irans-polyacetylene modeled by a dimerized 
tight binding chain. Interactions play no role in this con- 
struction and the many-body wavefunction indeed can 
be built as a Slater determinant of single-particle elec- 
tron states. In this case fractionalization is a many-body 
effect in that it is a cooperative phenomenon of a large 
number of electrons but it arises from the geometry of 
the system rather than strong correlations. 

In three dimensions, Jackiw and Rebbi Q argued that 
it is consistent to assign a half-integer fermion number 
to a monopole in the Yang-Mills gauge field coupled to 
relativistic Dirac fermions. In essence, this implies that a 
half-fermion is permanently bound to such a monopole. 
Again, interactions between fermions play no role in this 
construction and the half-fermion emerges from a many- 
body wavefunction that is a Slater determinant of single- 



particle states. It has been pointed out very recently [§| 
that a two-dimensional version of the physics underly- 
ing the above effect could be observable under certain 
conditions in graphcnc. 

An appealing feature of weakly interacting systems is 
that their properties can be often understood from es- 
sentially exact calculations, which is rarely possible for 
strongly correlated systems in dimension greater than 
one. It is thus desirable to study weakly interacting 
systems with exotic excitations if there is a chance that 
they can be realized in a laboratory. Our proposal for 
such a system is depicted in Fig. [TJa). It consists of a 
type-II superconducting film grown on top of a semicon- 
ductor heterostructure (such as GaAs/AlGaAs) hosting 
a 2D electron gas. When a magnetic field is applied per- 
pendicular to this planar device, an Abrikosov lattice of 
vortices forms in the superconductor. This results in a 
periodic modulation of the magnetic field in the 2DEG 
with one-half of a magnetic flux quantum <f>o = hc/e per 
unit cell. Devices like this have in fact been fabricated 
previously [ll|, [H, and studied theoretically [3, [H| . 

A key new aspect of our proposal is a periodic array of 
pinning sites imprinted on the superconductor. Pinning 
sites are regions of weakened superconductivity, which 
attract vortex cores and can be created in a variety of 
ways (l6| . When the number of vortices in the sample 
equals the number of pinning sites, each pin is occupied 
by a vortex and a periodic vortex lattice arises. The 
Magnetic field strength at which this happens is called 
the matching field Bm- When the field is slightly below 
(above) Bm, vacancies (interstitials) appear in the vortex 
lattice as illustrated in Fig. [T^b). Such defects in the 
vortex lattice then lead to a localized deficit or surplus 
of magnetic flux quantized in the units of $o/2 shown in 

Fig. nrc). 

Our principal claim, which we justify below by a gen- 
eral argument and a detailed computation, is this: When 
the lowest Landau level of the 2DEG is filled by spin- 
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FIG. 1: (a) The schematic of the proposed device. Side view 
of the superconductor-2DEG heterostructure in an applied 
magnetic field B. (b) Top view. Open circles represent pin- 
ning sites in the superconductor while crosses denote vortices. 
Throughout this paper we consider, for simplicity, the square 
vortex lattice. All of our conclusions, however, remain valid 
for any arbitrary periodic Bravais lattice, including the most 
natural triangular case, (c) The magnetic field produced by 
a periodic square lattice of vortices with one vacancy (lower 
left corner) and one interstitial (upper right corner). The field 
distribution is calculated from a simple London model fiol ] 
with the penetration depth A = a and the vortex core cutoff 
£ = 0.05A. (d) Thought experiment: creating an interstitial 
by adiabatically threading flux through the system. 



polarized electrons, i.e. at filling fraction v = 1, a defect 
in the vortex lattice binds fractional electric charge in the 
2DEG with the exact value — e/2 for a vacancy and +e/2 
for an interstitial. Such a defect constitutes a literal real- 
ization of Wilczek's model of an anyon as a bound state 
of magnetic flux and electric charge [TtJ • The exchange 
phase of this anyon 7r/4. 

The flux deficit/surplus can be, at least in principle, 
created, detected, and manipulated by the suite of scan- 
ning SQUID and Hall bar probes that have been devel- 
oped over the years to image [H, [ljj [2(| and manipulate 
[2lL |22| vortices in superconducting films and crystals. 
The proposed device could thus offer a unique opportu- 
nity to manipulate anyons with an unprecedented degree 
of control and advance the quest for a topologically pro- 
tected quantum computer. 

In order to execute a useful quantum computation, 
non-Abelian anyons are needed We remark that 

the device pictured in Fig. Q] could prove useful for the 
manipulation of such non-Abelian anyons if the 2DEG 
were put into the v = | FQH state, which is be- 
lieved to be described by the Moore-Read "Pfaffian" state 
[IS HH [HI- In that situation, vacancies and intersti- 
tials should bind non-Abelian anyons and the scanning 
SQUID probe could be used to locate these and possibly 
perform the braiding operations required to implement a 
fault tolerant quantum algorithm. 

II. ANYONS BOUND TO VACANCIES AND 
INTERSTITIALS: GENERAL ARGUMENTS 

So why does a vacancy bind fractional charge? The 
answer can be given at several levels of sophistication 
and is essentially implicit in the body of work on the 
FQH effect. The simplest way to see that this is true 
is to recall that electrons in the quantum Hall state at 
v = 1 form an incompressible fluid. This means that in 
the ground state the electron charge density p(r) tracks 
the magnetic field strength B(r) in accordance with the 
defining relation 



where N e and iV$ denote the total number of electrons 
and flux quanta in the sample, respectively. In an in- 
compressible state the above relation remains valid, to 
a good approximation, even for a spatially varying mag- 
netic field, as long as the variation is weak. In our device, 
a slowly varying field occurs in the limit when the mag- 
netic penetration depth A 3> a, a being the intervortex 
distance. In this limit, one obtains 

p ( T )~ ve W (2 ) 



This relation, which becomes exact when integrated over 
the region containing the vacancy, implies that at v = 1 
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a deficit of one half flux quantum produces a deficit of 
one half of an electron in the 2DEG. 

To see more clearly how the charge deficit /surplus 
comes about we may employ a version of Laughlin's argu- 
ment |26| similar to that used by Halperin to deduce the 
existence of edge states in integer quantum Hall liquids 
[2^ . The idea is to start from a perfect vortex lattice 
and imagine creating an interstitial by adiabatically tun- 
ing the flux of the extra vortex from zero to $o/2- We 
again consider the limit )>a. 

According to Faraday's law, VxE = — (l/c)(9B/9f), 
the time dependent magnetic field associated with in- 
troducing the extra vortex induces an electric field with 
concentric field lines as shown in Fig.QJd). In the quan- 
tum Hall liquid such an electric field produces current 
strictly perpendicular to the field, 



j = <T X y(E X Z) 



(3) 



where a xy = ve 2 jh is the quantized Hall conductance. 
Using the continuity equation for the electric charge, one 
can integrate this current to find the rate of change of 
the total charge Q in the region bounded by an arbitrary 
closed contour C enclosing the flux, 



dQ 
dt 



dl n • j = —a. 



rfl-E, 



(4) 



where n is a unit vector normal to the contour. Using 
Stokes theorem, the last integral can then be writen as 
an integral of V x E over an area bounded by C. This 
we can evaluate with help of Faraday's law to obtain 
dQ/dt = (a xy /c)d&/dt. It follows that the charge 5Q 
accumulated inside C during the process of adiabatic flux 
insertion is 



SQ = ev 



<5$ 



(5) 



For v — 1 and <5<I> = ±<&o/2 the accumulated charge is 
±e/2. It is to be noted that this argument depends only 
on the fundamental property ^ and the above result ^ 
should thus be valid unconditionally as long as the 2DEG 
remains in the quantum Hall plateau. 

Vacancies and interstitials, being bound states of 
charge and flux, (— e/2, — $ /2) and (e/2, <&o/2) respec- 
tively, acquire nontrivial Aharonov-Bohm phases upon 
adiabatic exchange, and can thus be thought of as 
anyons. The standard textbook counting procedure [28| . 
according to which a charge q taken around a flux $ ac- 
quires a phase 27r(g/e) (<&/$()), would imply a "semionic" 
exchange phase 9 = 7r/2 for both vacancies and intersti- 
tials, and a mutual exchange phase 9 m = — 7r/2. This 
conclusion, however, is incorrect for the following reason. 
The charge ±e/2 is built up from fermionic degrees of 
freedom and thus has attached to it half of an electron 
spin, or more precisely, half of the electron statistical 
phase, which needs to be taken into account when com- 
puting the exchange phase of the composite object. The 
simplest way to arrive at the correct result is to consider 



fusion rules for anyons. The relevant rule [28| states that 
the exchange phase O of a particle formed by combining 
n identical anyons with exchange phase 9 is O = n 2 9. 

Consider bringing together two interstitials. The re- 
sulting object is (e, $0)1 a bound state of an electron and 
a full flux quantum. The exchange phase of this object 
is 7r + 2n — 37r, the first term reflecting the intrinsic elec- 
tron exchange phase and the second the Aharonov-Bohm 
phase. But a 37r phase is equivalent to 7r, i.e. the bound 
state is a fermion. The above fusion rule with n — 2 thus 
determines the exchange phase of an interstitial to be 
9 = 7r/4. A similar argument leads to the identical result 
for a vacancy. These results will be confirmed below by 
explicit model calculations. 



III. CONTINUUM MODEL 

We first consider a very simple continuum model, 
which is exactly soluble and captures the essence of 
the effect discussed above. Consider the following 2- 
dimensional electron Hamiltonian, 



H = 



— (p--aY 

2m e V c / 



(0) 



with m e the electron mass, p the momentum operator in 
the x-y plane and A = Ao + <5A. Here Ao = \Ba{v x z) 
represents a uniform field in the z direction and 



5A = ^-4(r x z). 



(7) 



The total field seen by an electron is then 

B(r) = VxA = zB + zrj^ 5(r). 

The ^-function serves as a crude representation of the 
flux added by an interstitial (77 = 5) or removed by a 
vacancy (rj = — g) in the vortex lattice. 

It is straightforward to find the single particle eigen- 
states ipkm (r) of Ti. for arbitrary real rj by working in the 
polar coordinate basis. They are labeled by the principal 
quantum number k = 0,1,2,... and an integer angular 
momentum m. The allowed energy levels read 



-huj c [2k + 1 + \m — rj\ — (to — 77)] 



(8) 



where lo c = eBo/m e c is the cyclotron frequency. The 
spectrum ([5]) is depicted in Fig.[^a). For non-integer 77, 
in addition to the usual Landau levels, there also exist 
bound states associated with the extra flux. 

The eigenstates ipom in the lowest Landau level have a 
simple form 



iP 0m {z) = A m \z\-^z m e-^ 2 !\ 



(9) 



where z — {x+iy)/£s is the dimensionless complex coor- 
dinate, A m = [^|27r2 m - , 'r(l + m-77)]" 1 / 2 is the normal- 
ization constant and £b = y/hc/eB the magnetic length. 
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FIG. 2: (a) The energy levels (in units of hco c ) of Hamiltonian (6) as a function of fractional flux r\. The numerals denote 
degeneracy of bound state levels represented by dashed lines. Degeneracy of the extended states (solid lines) is proportional to 
the area of the system, (b) Charge density p(r), in units of e/2ir£%, for TV = 100 as a function of distance r from the origin. 
Solid line represents density po(r) for uniform magnetic field, dashed line has a flux tube carrying — <I>o/2 at the origin, (c) The 
excess charge density 5p(r) = po{r) — p(r). Inset shows the accumulated excess charge SQ(r) = 2n JJ r'dr'8p(r') in units of e. 



If we fill the lowest Landau level with electrons, then the 
many-body wavefunction can be constructed as a Slater 
determinant of ip0m{ z i)i where Z{ is a complex coordinate 
of the i-th electron. The usual simplifications using the 
properties of Vandermonde determinants then give 

= i*r n& - *i) e- Ei lzi|2/4 - ( io ) 

i i<j 

We are interested in the charge density implied by the 
above many-body wavefunction, which is given as the ex- 
pectation value p(r) = (^ol/ol^o) of the electron density 
operator p = ej^i <5(r — r,-). For a droplet composed of 
N electrons occupying the N lowest angular momentum 
states we obtain 

JV-l 

p(r)= e ^|^ 0m (r)| 2 . (11) 

m=0 

The charge density obtained by numerically evaluating 
Eq. (jTTjl is displayed in Fig. [51 both for a uniform field 
{ji = 0) and the vacancy (r) = — i) placed at the origin. 
We observe that the inclusion of a vacancy leads to sup- 
pression of the charge density near the origin, the excess 
charge being deposited at the outer edge of the system. 
The integrated charge deficit SQ(r) quickly saturates at 
e/2, confirming the conclusion reached above by general 
argument. The vacancy in this model indeed binds half 
of an electron. 

From the many-body wavefunction (|10[) one can obtain 
the charge and the statistical angle associated with a va- 
cancy by the standard calculation of the Berry phase. 
To compute the Berry phase it is convenient to perform 
a gauge transformation 

SA^ 5A-^-X7p (12) 

where z = \z\e lv . The resulting vector potential can be 
seen to vanish everywhere, except for a cut along the 



positive x axis. The effect of the gauge transformation 
is to multiply the eigenstates (J5J) by e~ %r,ip which has the 
effect of replacing by z~ n . In this "string gauge" the 
entire phase information is encoded in the wavefunction. 
The many-body wavefunction for a vacancy (v = — |) 
located at the complex coordinate w can thus be written 
as 

= M« Hiz, - u0 1/a *o({*}), (13) 

i 

where * ({^}) = U t <M ~ Zj)e~^ |z ' |2/4 . 

If we now adiabatically carry the vacancy along a 
closed contour C the wavefunction acquires a purely ge- 
ometric phase 

1 (C)=i£dw(* w \-^* w ). (14) 

Calculation of 7(C) proceeds exactly as in Ref. [|J and 
yields 

where $ is the total magnetic flux enclosed by the con- 
tour C. This result coincides with the Aharonov-Bohm 
phase acquired by a charge —e/2 carried along a trajec- 
tory enclosing flux $, confirming once again our earlier 
result for the electric charge of the vacancy. 

For two vacancies, located at w a and Wb the many- 
body wavefunction reads 

^w a w b =Nw aWb Y[(Zi - W a ) 1/2 { Zl - W b ) 1/2 ^ . (16) 

i 

If we now carry w a along a contour C that encloses Wf,, 
the latter held static, then Eq. (fT4"]) gives the Berry phase 

^--"(sH)- (17) 
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The first term reflects the Aharonov-Bohm phase as be- 
fore. Since the encircling operation can be thought of 
as two consecutive exchanges we must interpret the sec- 
ond term as twice the exchange phase of the vacancy, 
6> = 7r/4. 

IV. LATTICE MODEL 

Thus far we have ignored the periodic variation of the 
magnetic field induced by the vortex lattice. A question 
thus arises whether this could affect the emergence of 
anyons in realistic systems. To see that this is not the 
case we now consider a model which describes the limit 
of strong field variation and also includes the effect of a 
Zeeman interaction. 

If the effective Zeeman coupling in a 2DEG is suffi- 
ciently strong, then, in addition to a quasi-periodic vector 
potential, the electrons also feel a periodic scalar poten- 
tial. This effect becomes important in diluted magnetic 
semiconductors, such as Gai-^Mn^As or Cdi-^Mn^Te, 
where the effective gyro-magnetic ratio can be of the or- 
der ~ 10 2 - 10 3 . As argued in Ref. [HI, in the limit 
A <C a, electrons in a 2DEG become almost localized 
near the maxima of the Zeeman field and can be de- 
scribed by a tight binding model with one-half magnetic 
flux quantum per unit cell. 

This motivates our consideration of a tight binding 
Hamiltonian 

H = - )j /,,<'" <y, + h.c.) + J2 e '* c K (18) 

ij i 

where cj creates an electron on a site of a square lat- 
tice, and Uj and e, are real hopping amplitudes and on- 
site energies, respectively. The effect of magnetic field is 
included by the Peierls phase factors 

% = !f A ' d1 ' (W) 

In the following we consider the case of nearest and next- 
nearest neighbor hopping amplitudes t and t' respec- 
tively, but note the addition of longer range hoppings 
does not affect our results. 

In the absence of defects we have q = e =const and the 
magnetic flux through each plaquette is $o/2- The spec- 
trum of the Hamiltonian (fT8|) consists of two branches, 
reflecting the fact that the unit cell now contains two 
lattice sites. In the Landau gauge, A = B$yx with 
-Bo = $o/2a 2 , the allowed energies read 

1 /2 

i?k = e ± 2t [cos 2 k x + cos 2 k y + i'y' 2 sin 2 k x sin 2 k y ] , 

(20) 

where k is a wavevector drawn from the first magnetic 
Brillouin zone, 7 = t' /t, and we have set a = 1. For t' = 
the spectrum exhibits two Dirac points at (±tt/2, 7r/2). 
Inclusion of next-nearest neighbor hopping breaks the 
time reversal symmetry and produces a gap at these two 



points, as illustrated in Fig. [3ja). We are interested in 
the half-filled case with the lower band filled and the up- 
per band empty. This situation is analogous to the v = 1 
continuum case with Hall conductance a xy — e 2 /h dis- 
cussed previously. 

A vacancy in the vortex lattice removes flux $o/2 from 
the plaquettes adjacent to the affected site, modifies its 
on-site energy q and possibly also those of the nearby 
sites. While the total flux removed is exactly < I > o/2, the 
change of the on-site energies dei depends on the details 
of the underlying microscopic model for the Zeeman in- 
teraction. We find that, as expected by the general argu- 
ment, the total charge deficit associated with a vacancy 
depends critically on the total flux removed but is inde- 
pendent of 5ti or the details of the flux distribution. 

We have solved the model specified by Hamiltonian 
(fTSf by exact numerical diagonalization for lattice sizes 
up to 50 x 50 sites. We have considered various configu- 
rations of 5ti. In all cases we found that the total charge 
induced by an interstitial/ vacancy is ±e/2, to within ma- 
chine accuracy. This result continues to hold even in the 
presence of randomness in the on-site energies, as long 
as the randomness is weak compared to the gap. This is 
illustrated in Fig. [3jb-f ) . 



V. SUMMARY AND CONCLUSIONS 

Our foregoing discussion shows how two simple and 
well understood systems, a type-II superconductor and a 
2-dimcnsional electron gas in the integer quantum Hall 
regime, can exhibit some rather unusual properties when 
brought into close proximity. These include nontrivial ex- 
change statistics and charge fractionalization, properties 
that are normally thought of as hallmarks of strong corre- 
lation physics. In our proposed device these phenomena 
occur in what must be characterized as a weakly corre- 
lated system: the many-body wavefunctions of electrons 
in both subsystems are conventional Slater determinants 
of single-particle states. 

There are two key ingredients that give rise to the 
above unusual properties: (i) precise flux quantization 
in a type-II superconductor and, (ii) incompressibility of 
the 2DEG at integer filling v. The device displayed in 
Fig. Q] is designed to project a precisely quantized deficit 
or surplus of magnetic flux onto the 2DEG, which re- 
sponds by adjusting its many-body wavefunction to pro- 
duce a localized deficit or surplus of electric charge with 
precise value e/2. This phenomenon can be regarded as a 
consequence of gauge invariance and is therefore robust 
against any weak but otherwise arbitrary perturbation 
imposed on the system. 

A question that we must ask is to what extent are the 
fractional charges associated with the vacancies and in- 
terstitials legitimate quasiparticle excitations of the sys- 
tem. Surely, from the point of view of a 2DEG at integer 
filling, such objects are not natural excitations. Indeed, 
in our device they occur in response to the local change 
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FIG. 3: a) Band structure of the tight binding model on a square lattice with one half flux quantum threading each plaquette 
and t' = 0.05t. b,c) Charge density p at half filling with random on-site energies (Gaussian distribution with Ae = 0.05t') and 
a full flux quantum modeling the interstitial piercing the central plaquette. d) Charge density po with the same realization of 
randomness but without the extra flux. e,f) The charge density Sp — p — po induced by the interstitial. The induced charge 
SQ integrates to e/2 to within machine accuracy. 



in the magnetic flux, which is external to the 2DEG. The 
key point is that fractional charges bound to defects are 
excitations of the entire system in the same sense as frac- 
tional charges bound to domain walls are excitations of 
a dimerized polyacetylene chain 0]. 

Fractional charges should be relatively easy to detect 
in a transport experiment in the regime where the longi- 
tudinal conductivity of the 2DEG is furnished solely by 
the e/2 defects. Fractional statistics, on the other hand, 
will be more difficult to establish since this requires an 
interference experiment; in essence we need a defect to 
quantum-delocalizc and interfere with itself after encir- 
cling another defect. As noted recently, non-trivial statis- 
tics might be easier to detect in the v — | non-Abelian 
case |29l l30j. 

Aside from demonstrating the possibility of uncon- 
ventional phenomena in a weakly interacting system, 
our proposed device could allow for the manipulation of 
anyons in the context of quantum information processing. 
Indeed, some schemes for fault tolerant topological quan- 



tum computation [J, Q involve execution of braiding 
and fusion operations on anyons. In the present setup 
these could be facilitated by the fact that our anyons 
are permanently bound to physical magnetic fluxes and 
well developed techniques exist to detect and manipulate 
these. As already remarked above if the 2DEG in our 
device were tuned into the v = | Moore-Read Pfaffian 
state then a defect in the vortex lattice would bind a 
quasiparticle with non-Abelian exchange statistics. Hav- 
ing a "handle" on such a quasiparticle in the form of the 
attached magnetic flux could prove useful for realizing 
the requisite braiding and fusion operations. 
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